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Is there a redshift cutoff for submillimetre galaxies? 
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ABSTRACT 

We present new optical and infrared photometry for a statistically complete sample of 
seven 1.1 mm selected sources with accurate coordinates. We determine photometric 
redshifts for four of the seven sources of 4.64, 4.54, 1.49 and 0.18. Of the other three 
sources two are undetected at optical wavelengths down to the limits of very deep 
Subaru and Canada-France-Hawaii Telescope images (~27 mag AB, i band) and the 
remaining source is obscured by a bright nearby galaxy. The sources with the highest 
redshifts are at higher redshifts than all but one of the ^200 sources taken from the 
largest recent 850 ^m surveys, which may indicate that 1.1 mm surveys are more 
efficient at finding sources at very high redshifts than 850 fj,m surveys. 

We investigate the evolution of the number density with redshift of our sample 
using a banded Ve/Va analysis and find no evidence for a redshift cutoff, although the 
number of sources is very small. We also perform the same analysis on a statistically 
complete sample of 38 galaxies selected at 850/^m from the GOODS-N field and find 
evidence for a drop-off in the number density beyond z 1 and 2, confirming the earlier 
conclusion of Wall, Pope & Scott (2008). We also find evidence for the existence of 
two differently evolving sub-populations separated in luminosity, with the drop-off in 
density for the low- luminosity sources occurring at a lower redshift. 

Key words: galaxies:distances and redshifts, galaxies:evolution, galaxies:high- 
redshift, submillimetre, infrared: galaxies, galaxies:statistics 



1 INTRODUCTION 



Submillimetre (submm) galaxies (SMGs), first detected at 
850 nm with the Su bmillimetre Common User Bolome- 
ter Array (SCUBA) jHolland et al. I Il999 l'). are a signifi- 
cant population of high redshift star f orming galaxies (e.g., 
iHughes et al. ll998l : [Blain et al. |2003 ). They are believed to 
be dust enshrouded galaxies undergoin g prodi gious levels o f 
star formation (e.g.. iHughes et al. Ill9 98: Ealcs et al. Ill999l ) 
in which the optical/UV radiation emitted by the stars is 
absorbed by the dust and then re-emitted in the submm. 
Star formation rates in excess of 1000 Meyr"^ have been 
inferred (|Scott et al. Il2003 ). much higher than locally. The 
galaxies in these samples have been found to account for 
up to one tenth of the total far-i nfrared/submm energetic 
background (e.g.. lOve et al. II2007I ) and many authors have 
argued that these galaxies are the progenitor s for the el- 
liptical galaxies we see i n the local Unive rse (Eales et al. I 
[l999; Scott ot al. 2003; iDunne. Eales fc E dmunds 2003). 
Thus understanding the nature of these sources is of great 
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importance for the understanding of galaxy formation and 
evolution as a whole. 

Observations of SMGs at ~1 mm benefit from a neg- 
ative K-correction out to high redshifts due the shape of 
their spectral energy distribution (SED). As the redshift of 
an SMC increases, the peak of its rest-frame SED moves to- 
ward the observed waveband, offsetting the dimming caused 
by the increasing luminosity distance. This fact accounts for 
the surprising ability of SCUBA to find large numbers of 
high-redshift galaxies. 

The large amount of dust responsible for the strong 
submm emission gives rise to high levels of attenuation in the 
optical. This in conjunction with the poor angular resolution 
of single dish submm facilities makes the cross identification 
of SMGs at different wavelengths difficult. Moreover, even 
when an optical counterpart can be identified, the high levels 
of dust attenuation makes the determination of a spectro- 
scopic redshift difficult. As such we are currently unable to 
determine spectroscopic redshifts for the majority of SMGs. 
The strong correlation between dust emission and radio 
emission which appears to hold true in both the low-redshif t 
and high-redshift universe (|Vlahakis. Eales fc Dunne Il2007h 
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has been useful for both identifying the counterparts and 
estimating redshifts. Due to the low surface density of ra- 
dio sources on the sky, the probability of the radio counter- 
part being coincidental with the submm source by chance 
is small. Due to the high positional accuracy of radio obser- 
vations, it is then possible to identify the optical counter- 
part and retrieve a spectroscopic redshift. It is also possible 
to estimat e the redshift using t he ratio of radio to submm 
flux (e.g.. iHuehes et al.1 Il998l : ICarilh fc YurTI Il999l . I2OO0I : 
ISma,il et al. II2OO0I). 

Chapman et al. I l|2005l '). u sing the Low R esolution 
Imaging Spectrograph (LRIS) (jOke et al. I 1 19951 ) on the 
Keck I telescope, managed to obtain spectroscopic redshifts 
for a total of 73 radio-identified SMGs with a median 850 
^m flux of 5.3 mjy. The galaxies in this sample were found 
to lie at a median redshift oi z = 2.2 out to a maximum value 
of Zmax = 3.6. However, the K-correction which allows us to 
detect high-redshift SMGs does not similarly benefit their 
radio fluxes and so radio identified SMGs are subjected to 
a ra dio selection effect which limits redshifts to 2 < 3. 

IPope et al. I |2006l). produced the first complete (i.e. not 
requiring radio IDs) sample of 850 fim selected SMGs that 
has close to 100% redshifts. The sample consists of 35 galax- 
ies, 21 with secure optical counterparts and 12 with tentative 
optical counterparts, and its completeness means that unlike 
previous surveys it is not biased towards low-z sources. The 
median redshift determined for this s a mple is z ~ 2.2. Us- 
ing this sample, IWall. Pope fc Scott I (|2008l ) examined the 
epoch dependency of the number density of SMGs. They 
found an apparent redshift cutoff a.t z > 3 with further ev- 
idence for two separately evolving populations, divided by 
luminosity. However this result was based on calculations 
using a single model galaxy SED. Since the predicted rela- 
tionship between submm flux-density and redshift depends 
strongly on the assumed SED, one of the aims of this paper 
is to re-examine their conclusion using a range of empirical 
SEDs. 

There have been a number of explanations for the lack 
of high-redshift SMGs. Given that dust is thought to form 
in the atmospheres of highly evolved stars, it is possible that 
at high r edshifts simply not enough t ime has passed for dust 
to form (|Morgan fc Edmunds II2003I ). Observations of high- 
redsh ift quasars have however det ected high levels of dust 
(e.g., IPriddev fc McMahoiTI 120011 ). suggesting that this is 
not the explanation. Another possible explanation is that 
ther e are fewer large st ar-forming galaxies at high redshifts. 

lEales et al. I l|2003l ) presented evidence that SMGs typ- 
ically have low values for the ratio of the 850 /im to 1200 
^m fluxes compared to that expected from a low redshift 
galaxy. One possible explanation is that these sources are 
at very high redshifts. If this is true, then observations at 
1.1 mm would be better at detecting SMGs at the highest 
redshifts than observations at 850 fim. A new complete sam- 
ple of 1.1 mm selecte d SMGs located in the COSMOS held 
( Scoville et aTlbOOTl ) has been compiled bv I Younger et slL] 
( 20071 ). The sources were selected init i ally at 1.1 mm with 
the AzTEC camera (|Scott et al. II2OO8I : IWilson et ani2008l ) 
on the JCMT. The resultant catalogue consists of 44 sources 
with S/N> 3.5(7, 10 of which are robust with S/N> 5a. 
Follow up observations bv I Younger et al. I (|2007h were then 
made with the Submillimetre Array (SMA) at 890 fim for 
the 7 highest significance AzTEC sources, allowing their po- 



sitions to be determined with an accuracy of ~ 0.2". The 
COSMOS field offers a wealth of data over a great num- 
ber of wavebands including the optical and infrared. Thus 
the high positional accuracy allows for the identification of 
optical counterparts and hence the determination of photo- 
metric redshifts. Of the seven AzTEC sources imaged with 
the SMA six have IRAC counterparts, and one source is ob- 
scured by a nearby bright galaxy. Using deep Hubble Space 
Telescope (HST) imagi ng acquired with th e Adva nced Cam- 
era for Surveys (ACS). iKoekemoer et al."1 l|2007l ) found op- 
tical counterpart candidates for only three of these sources. 

The main aim of this paper is to carry out a deeper 
search for the optical counterparts for the AzTEC sources. 
We give photometry from deep Subaru and Canada-France- 
Hawaii- Telescope (CFHT) imaging and find one new possi- 
ble optical ID. We estimate photometric redshifts for the 
AzTEC sources using the HyperZ photom etric redshift 
package (|Bolzonella. Miralles fc Pello [ I2OO0I ). Throughout 
this work we employ a concordance cosmological model with 
^itotai = 1, rim = 0.3, = 0.7 and Ho = 75 kms"^Mpc"\ 
All magnitudes quoted are AB. 



2 NEW IMAGES AND PHOTOMETRY FOR 
THE AZTEC SAMPLE 

We searched for optical counterparts and measured new pho- 
tometry using deep SubarrQ, CFHT and IRAC images of 
the AzTEC sources. The IRAC and Subaru images are the 
pubhc ly available COSMO S images taken by the COSMOS 
team l|Scoville et al. II2OO7I ). The CFHT images are taken 
from the CFHT Deep Legacy Survey. The images we used 
were taken using the CFHT gu, tm, im, zm, Subaru Bj, Vj, 
T+, i+, z+ and IRAC channel 1 and 2 filters to average 3a 
depths of approximately 28.4, 27.9, 27.6, 26.5, 29.0, 28.2, 
28.3, 27.7, 26.4, 24.1 and 23.6 mags respectively. 

We searched the i-band images (figure[2ll at the SMA co- 
ordinates. We find bright i-band counterparts for AzTEC 1, 

3 and 7, all of which were previously known. We also find 
a faint i-band counterpart for AzTEC5 at the SMA coordi- 
nates. We find no objects directly at the SMA coordinates 
for AzTEC2, but there is a bright object offset from this 
position by 3", meaning that the magnitude limits of this 
SMC are not useful. 

We find no optical counterparts directly at the SMA co- 
ordinates for AzTEC4 and 6 in the Subaru and CFHT imag- 
ing. For the latter source, however, there is a bright i-band 
counterpart offset from the SMA position by ~0.6" (~3cr) 
which could be AzTEC6's counterpart or the true counter- 
part may be too faint to see. There is also a faint i-band 
source, offset from AzTEC4's SMA position by ~0.8" (~4cr), 
in the Subaru imaging. For these two sources we added the 
i band CFHT and Subaru images, inversely weighting the 
images by the square of the noise, in order to try and detect 
any very faint possible counterparts. The stacked i-band im- 
ages for AzTEC4 and 6 are shown in figure [5] We still do not 
find counterparts at the SMA positions for AzTEC4 and 6 



^ An additional uncertainty of 0.3 mags in the Subaru Bj band 
magnitudes is taken into account in this photometry due to the 
possibiUty of a red leak or a shift in the blue cutoff of this filter. 
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and given the good coincidence between the SMA and op- 
tical positions for the other AzTEC sources we tentatively 
conclude that the true counterparts have not yet been de- 
tected. 

The typical full-width half-maximum (FWHM) of the 
optical and IRAC channel 1 and 2 point spread functions 
(PSFs) are ~0.8", 1.66" and 1.72" respectively. Magnitudes 
were determined manually by placing apertures onto the im- 
ages, ensuring that the aperture was large enough to contain 
as much of the emission from the galaxy as possible with- 
out including any emission from neighboring objects. Thus 
the sizes of the optical apertures vary from source to source, 
although are constant for a given source. We used larger 
apertures for the IRAC sources due to the images having a 
larger PSF, but use the procedure outlined below to correct 
for this. 

Due to the difference in the PSF between the optical 
and IRAC images as well as the difference in the aperture 
sizes, a small correction needed to be applied to the IRAC 
magnitudes before they could be used in conjunction with 
optical magnitudes to determine a photometric redshift. We 
corrected IRAC magnitudes by firstly fitting a 2D Gaussian 
to the IRAC source. We then scaled it to have the FWHM 
it would have had if observed with CFHT/Subaru. The fiux 
was then computed using the scaled Gaussian and new aper- 
ture size. All the corrections applied in this work increase 
the IRAC magnitudes, and the more extended the source 
the greater the correction. Corrections range from 0.01 to 
0.46 magnitudes. 

The new photometry is summarized in table [1] and the 
details of the individual objects are discussed below. 

2.1 Notes on Individual Objects 

AzTECl- J095942. 86+022938. 2- AzTECl is the 
brightest submm source in the sample with fiuxes of 
-FsgoAim ~ 15.6 ± 1.1 mjy and F\.imm = 10.7 ± 1.3 mjy. 
There is a bright i-band object located directly at the SMA 
position. Optical fluxes are measured using an aperture 
1.94" in diameter and the source is detected in the Subaru 
i-l- band at 25.11±0.03 mag whic h is in agreement wit h 
the HST i-band magnitude given in lYounger et al~l l|2007h . 
There is some disagreement within the same wavebands 
between the CFHT and Subaru and photometry (table [T]), 
but the discrepancy is small enough that the photometric 
redshift is not significantly affected. IRAC magnitudes were 
measured using an aperture of diameter 4.45". Only a small 
correction was applied to the IRAC magnitudes: -1-0.01 mag 
in both IRAC channel 1 and 2. 

AzTEC2- J100008. 05+022612.2- AzTEC2 is detected 
in the submm with fluxes -Fggoum = (12.4 ± 1.0) mJy and 
Fi.imm = (9.0 ± 1.3) mJy. No objects are found directly 
at the SMA coordinates, but there is a bright object offset 
from SMA position by 3" . Thus the limit on the magnitude 
of the optical counterpart is not very useful. 

AzTECS- J100020. 70+023520. 5- AzTEC3 is detected 
in the submm with fluxes -FsgoMm ~ 8.7 ± 1.5 mJy and 
Fi.imm = 7.6 ± 1.2 mJy. There is a bright i-band object 
located at the SMA coordinates as well as three companion 
objects offset by between 1" and 2". Since SMGs often 



seem to consist of multiple components (jlvison et al. |[l998h 
it is possible that these companion objects are also part of 
AzTEC3. However, since six of the seven AzTEC sources 
are detected in the IRAC bands, it seems likely that if 
the companion objects are part of the same galaxy then 
they should also be contributing to the IRAC emission. 
We attempted to determine whether this is the case by 
convolving the Subaru image with the IRAC beam and 
comparing the FWHM of the IRAC source with that of the 
convolved Subaru image. We find that the FWHM of the 
convolved image is ~4.4". The FWHM of the IRAC 3.6 
Hm image is ~2.87", suggesting that the 3.6 /xm emission 
is associated only with the central object. Optical fiuxes 
were measured using a aperture of diameter 1.26" and the 
source is detected in the Subaru i-f band at 26.18±0.08 
mag which is in agreement with the HST i-band magnitude. 
The CFHT and Subaru magnitudes within the same bands 
are consistent with each other. IRAC magnitudes were 
measured using an aperture of diameter 4.80". A correction 
of -1-0.2 mag was applied to the IRAC magnitudes in 
channels 1 and 2. 

AzTEC4- J095931. 72+023044.0- AzTEC4 is detected 
in the submm with fiuxes Fs90iJ.m = 14.4 ±1.9 mJy and 
Fi.imm = 6.8 ± 1.3 mJy. We find a tentative i-band 
counterpart, offset from the SMA position by 0.8" (~3(t), 
in the Subaru image with a magnitude of 27.43±0.13. 
In the combined image (see above), the counterpart can 
be seen more clearly and has a magnitude of 26.99±0.18 
in a 2.57" diameter aperture. However we find it is too 
faint to detect in our other Subaru and CFHT bands. 
IRAC magnitudes were measured using an aperture of 
diameter 4.80". Corrections of +0.11 and -f0.04 mags 
were applied to the IRAC magnitudes in channels 1 and 2 
respectively. Because of the good agreement between the 
SMA the optical positions for the other AzTEC sources 
we tentatively conclude that this is not the true counterpart. 

AzTECS- J100019. 75+023204. 4- AzTEC5 is detected 
in the submm with fiuxe s Js90um = 9.3 ± 1.3 mJy and 
Fl,lmm — 7.6 ± 1.3 mJy. lYounger et aTl (|2007D found no 
optical counterpart in ACS imaging, but we find a faint 
Subaru source at the SMA coordinates with a Subaru i+ 
band magnitude of 26.74±0.13, measured in an aperture of 
diameter 1.68". The CFHT and Subaru magnitudes within 
the same bands are consistent with each other. IRAC 
magnitudes were measured using a aperture of diameter 
4.80". Corrections of -1-0.46 and -1-0.14 mag were applied to 
the IRAC magnitudes in channels 1 and 2 respectively. 

AzTECe- J100006. 50+023837.7- AzTEC6 is detected 
in the submm with fiuxes Js90um = 8.6 ± 1.3 mJy and 
Fi.imm = 7.9 ± 1.2 mJy. lYounger et al~l l|2007l ) find no 
optical counterpart in ACS imaging. In CFHT and Subaru 
imaging we find no source directly at the SMA coordinates, 
but we do find a source offset from the SMA position 
by ~0.6" (~3(j). This could therefore be the optical 
counterpart, or the true counterpart may be too faint to 
detect. The source offset from the SMA position has a 
Subaru i-|- magnitude of 25.38±0.04 magnitudes, measured 
in an aperture of diameter 1.59". The CFHT and Subaru 
magnitudes within the same bands are consistent with each 
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Figure 1. Subaru i+ band cutouts for AzTECl to 7, with the exception of AzTEC4 and 6, which are the combined CFHT and Subaru 
i-band images. Each image has a field of view of 15.3"xl5.3" and a scale of 0.15"/pixel. The SMA coordinates of each source are 
highlighted by cross-hairs and the optical counterpart (including the objects offset from AzTEC4' and 6s SMA coordinates by ~0.8" and 
~0.6" respectively, see text) is circled with an aperture the size of which was used in its photometry. SMA coordinates are accurate to 
~0.2". 



other. IRAC magnitudes were measured using an aperture 
of diameter 5.88". A correction of +0.13 mag is applied 
to the IRAC magnitudes in channels 1 and 2. Because of 
the good agreement between the SMA the optical positions 
for the other AzTEC sources we tentatively conclude that 
this is not the true counterpart, although we do estimate a 
photometric redshift for it. 

AzTEC7- J100018.06+024830.5- AzTECT is detected 
in the submm with fluxes Fsoojim = 12.0 ± 1.5 mjy and 
Fi.imm = 8.3 ± 1.4 m,Jy. We find an optical counterpart 
with a disturbed morphology at the SMA coordinates which 
could be a system of merging galaxies. Optical fluxes were 



measured by placing an aperture of diameter 2.87" over the 
whole of the system. The source is detected in the Subaru i+ 
band at 24.20±0.04 mag. IRAC magnitudes were measured 
using an aperture of diameter 6.12". Corrections of +0.08 
and +0.02 mag were applied to the IRAC magnitudes in 
channels 1 and 2 respectively. 
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Figure 1 - continued 



3 ESTIMATED REDSHIFTS FOR THE AZTEC 
SAMPLE 

Photometric redshifts were determined by apply- 
ing the photometric re d shift package, HyperZ 
l|BoIzonella. Miralles fc Pello I |2000| '). to om- 11 band 
photometry (Subaru:B, V, r+, i+, z+; CFHT: gyi, m, 
«M, zm\ IRAC: 3.6 /xm, 4.5 ^m). The spectra used for 
fitting in this w ork are taken from the set compiled by 
iDve et al~l ()2008l 'l. which is optimized for the determination 
of photometric r edshifts wh e n inc luding filters in the 
near/mid-infrared. iDve et al. I (|2008l ) compared the photo- 
metric redshifts determined using these spectral templates 
with those determined using synthetic spectra constructed 
from the best-fit star formation history for their sample of 
60 SCUBA sources. Since these methods are completely 
independent and the redshifts found using both sets of 



templates were found to be in good agreement, we assume 
that our template set is adequate. 

We varied the redshift in t he range z = t o 10. W e em- 
ployed the reddening regime of ICalzetti et al. ] (|2000l ). with 
Av allowed to vary in the range Ay = to 5 in steps of 
0.2. We used a minimum photometric error of 0.05 magni- 
tudes for each band. For wavebands in which we have no 
detection we took the fiux of the source to be zero with 
a 1(7 error equal to the sensitivity of the detector in that 
waveband. The photometric redshifts obtained are listed in 
table H 

The median redshift of the sample is 2.7 which is some- 
what high er than the media n redsh ift, 2.2, of the sample pre- 
sented by IChapman et al. 1 l|2005h . The maximum redshift 
found is 4.64 and the minimum redshift found is 0.18. Com- 
paring the redshift dis tribution of this sample t o that of the 
samples presented in I Chapman et al. ] l|2005l ). iPope et al. I 
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AzTECl 


AzTEC2 


AzTEC3 


AzTEC4 


AzTEC5 


AzTECe 


AzTEC 7 


RA 


nQ-5q-42 86 


10 00 08 05 


10-00-20 70 


09-59-31 72 


10-00-19 75 


10-00 06 50 


10 00 18 06 


Dec 




_i_n2-2fi-1 2 2 


-I-02 S5-20 5 


-I-02-30-44 


-I-02-32 04 4 


4-02-38-37 7 


+02-48-30 5 




1.94" 




1.26" 


2.57" 


1.68" 


...yi.Oi) ) 


2.87" 


vn^ 


>28.88 




>29.14 


>28.67 


2« 80+0 47 


^28 Q8('25 80+0 30"! 




my 


27.13±0.17 




28.77±0.75 


>28.21 


28.76±0.62 


>28. 32(25.67+0.04) 


25.10+0.06 


mr+ 


26.21±0.06 




27.39±0.22 


>28.02 


27.07±0.13 


>28. 47(25.77+0.03) 


24.97+0.05 


mi+ 


25.11±0.03 




26.18±0.08 


(27.43±0.13) 


26.74±0.13 


>27.97(25.38+0.04) 


24.20+0.04 




25.02±0.02 




25.58±0.15 


>26.50 


26.07±0.22 


>26. 79(24.80+0.06) 


23.65+0.07 


9M 


>28.12 




>28.71 


>27.93 


>28.41 


>28. 42(26.16+0.05) 


N/A 


rm 


26.54±0.15 




27.13±0.24 


>27.46 


27.15±0.20 


>27.96(25.61+0.05) 


N/A 


IM 


25.25±0.05 




26.30±0.12 


>26.19 


26.50±0.13 


>27.81(25.42+0.05) 


N/A 




25.11±0.13 




25.69±0.22 


>26.27 


26.46±0.38 


>26. 06(24.87+0.08) 


N/A 


IRAC Ap. Size 


4.45" 




4.80" 


4.80" 


4.80" 


5.88" 


6.12" 




23.40±0.07 




23.72±0.11 


22.11±0.04 


23.24±0.08 


24.13+0.25 


20.63+0.01 




23.08±0.08 




22.98±0.12 


22.15±0.04 


22.31±0.06 


23.50+0.27 


20.15+0.02 



Table 1. Photometry for the AzTEC sources, given in AB magnitudes. First two rows give the SMA co-ordinates. Aperture sizes are 
the diameters used for measuring optical and IRAC magnitudes. The IRAC magnitudes are corrected to take into account the difference 
in the seeing and aperture sizes for the IRAC and optical imaging (see text). No optical counterparts were found for AzTEC2. The only 
nearby optical counterparts for AzTEC4 and 6 are offset from their SMA positions by ^4(7 and ^3(7 respectively. We give the photometry 
for these objects in parentheses. 



ID 



Notes 



AzTECl 
AzTEC2 
AzTEC3 

AzTEC4 
AzTEC5 



4.64 + 0.06 



4.54 + 0.10 



1.49 + 0.10 



AzTEC7 0.18 + 0.01 



1.537 



2.196 



1.488 



AzTEC6 (2.09 + 0.01) (6.172) 



7.021 



No optical counterpart. 
There is a secondary chi-squared minimum at the lower 
redshift of 2 ~ 0.4 with a chi-squared fit value of ^^3.5. 
Nearest counterpart only detected in one optical band. 
There is a secondary chi-squared minimum at the higher 
redshift of z ~ 4 with a chi-squared fit value of ~4. 
The redshift and chi-squared values are for the optical source 
offset from AzTEC6's SMA position. The chi-squared fit to 
this source is much poorer compared to the others in the sample. 
This may further imply that the nearby optical counterpart we have 
selected is not the true counterpart to AzTEC6 and that the IRAC 

emission is unassociated with the optical emission. 
CFHT data not available. There are several other possible redshifts 
with chi-squared fit values of ~ 10 up to 2 ~ 2. Even the best 
chi-squared fit is still relatively poor however, which may be due 
to the unusual nature of the source. 



Table 2. The best photometric redshift fits for the sources with their minimum reduced x value, x 



Notes of interest on the 



photometric redshifts, including secondary fits, for each source are also given. We do not give the best fit SED type as typically for each 
source there are several SED types which fit equally well. Note that reduced chi-squared values given here are not those directly output 
by HyperZ, which takes the number of degrees of freedom as being the (number of filters— 1). This is true if only the redshift is allowed 
to vary. However, we are additionally allowing Av, SEDs type and the normalization to vary. Thus the correct number of degrees of 
freedom is given by (number of filters— 4). 



l|2006t ). iDve et al. I (|2008D and IClements et al. I (|2008D . we 
note that only one of the sources in this combined sample of 
~200 850 /im selected sources is at a comparably high red- 
shift as our two highest redshift sources, although this dif- 
ference is not significant when analyzed with a Kolmogorov- 
Smirnov test. However two of the other AzTEC sources are 
undetected to very faint limits in the i-band, and these facts 
may indicate that 1.1 mm surveys find more sources at very 
high redshifts than 850 /.tm surveys. 



4 A BANDED Ve/Va ANALYSIS 

IWall. Pope fc Scott I l|2008l ) examined a sample of 38 SMGs 



in the GOODS-N field and found evidence for a diminution 
in the space density of SMGs at redshifts 2; > 3. They also 
found evidence for two separately evolving sub-populations 
separated by luminosity. In this paper we present the results 
of our re-examination of this result using a banded Ve/Va 
analysis and a range of empirical SEDs rather than the the- 
oretical SED used by Wall et al. 

The most well known method of investigating the evo- 
lution of the space density of galaxie s with redshift is the 
{V/Vmaa/) test (Schmidt 1968; Rowa n-Robinson Ill968l ). V 
is the co-moving volume enclosed by the galaxy (that vol- 
ume which the field of view traces out in moving from a 
redshift of z = out to the galaxy) and Vmax is the vol- 
ume that would be enclosed by the galaxy were it pushed 
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Figure 2. The left hand column shows the photometric data points for the AzTEC sources with optical counterparts. The best spectral 
fits for the sources are overlaid. The right hand column shows the marginalized reduced distribution as a function of redshift. The 
AzTEC6 plots correspond to the nearby optically bright object. 



to the redshift at which its flux drops to the survey limit. 
This method encounters problems when a survey encloses 
two galaxy populations, one undergoing positive evolution, 
and the other negative. If we have a uniform distribution of 
galaxies in space, then we expect the value of {V /Vmax) to 
be 0.5±(12A'')"°-^, where iV is the number of sources in the 
sample. A value of (y/Vmax) > 0.5 then implies a concen- 
tration of sources toward the more distant regions of their 
accessible volume and a value of (V/Vmax) < 0.5 implies a 
deficit of sources at higher redshifts. Therefore if we have in 
our sample separate populations undergoing high levels of 
positive and negative evolution, then (V/Vmax) niay still be 
close to 0.5, incorrectly implying zero evolution. 



This pro blem can be solved by imp lementing instead a 
{Ve/Va) test (|Dunlop fc Peacock I ll990l ). This is effectively 
a banded version of the (V/Vmax) test. Ve, the effective vol- 
ume, is the volume enclosed between a minimum redshift 
ziom and the redshift of the galaxy. Va, the accessible vol- 
ume, is the volume enclosed between ziow and the redshift 
at which the galaxy's fiux drops below the sensitivity of the 
survey. By investigating the variation of (Ve/Va) with 2;;o„ 
we can distinguish between a positively evolving and a neg- 
atively evolving population. 

We investigated the evolution of the space density of 
the sample with redshift through the implementation of a 
(Ve/Va) test. Wall et al. based the k-correction necessary 
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Figure 2 

to calculate accessible volume on a single theoretical SED, 
whereas real galaxies have a range of SEDs. To investigate 
this, we carried out the (Ve/Va) analysis using two different 
assumptions about SEDs. We used the two ext reme two- 
component dust models of lDunne fc Ealesll|200lD . who pro- 
vided fits to the hottest and coldest local SMGs. The cold 
SED, based on NGC 958, contains dust at temperatures of 
20 and 44 K with a cold-to-hot dust mass ratio of 186:1. The 
hot SED, based on IR1525-I-36, contains dust at tempera- 
tures of 19 and 45 K with a cold-to-hot dust mass ratio of 
15:1. Figure|3l which shows the predicted flux versus redshift 
plot for the different models, shows the effect of using dif- 
ferent SED templates on the flux-redshift relation. The two 
SED types are normalized such that they produce a flux of 
1 mjy at a redshift of z = 1. 

We took the limiting flux of each source in the GOODS 
sample to be 3.5(j and measured {Ve/Va) for ziow = to 4 in 
steps of 0.1. We also separated sources into two samples of 
equal size according to luminosity. In doing this we are able 
to determine whether there are differences in the evolution 
of the two sub-populations. 

Our results for the 38 SMGs of Wall, Pope & Scott 
are shown in figure |4] We find good evidence for the exis- 
tence of a redshift cutoff at 2; > 1 for the hot SED, and 
slightly weaker evidence for a redshift cutoff at 2 > 2 for 
the cold SED. Dividing the sample into separate popula- 
tions of high and low luminosity sources shows differences 
in the evolution of the two populations. The low luminos- 
ity sources show much sharper redshift cutoffs whereas the 
evidence for redshift cutoffs in the high luminosity sources 
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Figure 3. Flux versus redshift for both the cold (solid line) and 
hot (dashed line) SEDs. Both SEDs are normalized such that they 
produce a flux of 1 mJy at a redshift of z = 1. 

is far more marginal . Thus we find evid ence to support the 
conclusions given in IWaU. Pope fc Sco tt (2008): there IS a 
redshift cutoff for the sample and that there is evidence for 
two separately evolving sub-populations 

An add i tiona l uncertainty about this results is that 
iPope et al. I (|2006| ) claim that only 60% of their identifica- 
tions are reliable. Therefore we also performed the (14/14) 
analysis only on sources with reliable identifications, the re- 
sults of which are shown in figure[S] Using these sources only, 
we still find good evidence for a redshift cutoff at z > 1 for 
the hot SED, and some marginal evidence for a cutoff at 
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Figure 4. The distribution of the values of {V^/Va) with zi^nj for the GOODS-N sample. Figures in the left hand column are for hot 
SEDs and figures in the right hand column are for cold SEDs. The sample is also separated into high and low luminosity sources. The 
dashed line denotes the position of (Ve/Va) = 0.5 on the graph, values above which imply a concentration of galaxies at higher redshifts 
and below which imply a concentration of galaxies at lower redshifts. The grey shaded region shows the area enclosed by the la error. 



z > 2 for the cold SED. However, we are unable to find 
any clear evidence for two separately evolving galaxy sub- 
populations, separated in luminosity, as the sample size is 
too small. 

However, by only taking into account the reliable identi- 
fications, we are probably biased towards optically brighter 
galaxies and therefore lower redshifts. We further investi- 



gated the effect of the unreliable identifications by putting 
four (roughly half) of the unreliable identifications ai z — 4 
and repeating the analysis (figure |6} . Doing this, we find 
that for hot SEDs our results are largely unaffected, with a 
relatively clear cutoff at redshifts higher than z = 1. How- 
ever for the cold SEDs we find that our results are strongly 
affected, with no clear redshift cutoff up to a redshift of z~3. 
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Figure 5. The distribution of the values of {Ve/Va) with ziow for the GOODS-N sample, using only the sources with reliable identifica- 
tions. Figures in the left hand column are for hot SEDs and figures in the right hand column are for cold SEDs. The dashed line denotes 
the position of {Ve/Va) = 0.5 on the graph, values above which imply a concentration of galaxies at higher redshifts and below which 
imply a concentration of galaxies at lower redshifts. The grey shaded region shows the area enclosed by the la error. 
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Figure 6. The distribution of the values of {Ve/Va) with z;q„ for the GOODS-N sample, where four (roughly half)of the unreliable 
identifications have been pushed to redshifts of 2 = 4. Figures in the left hand column are for hot SEDs and figures in the right hand 
column are for cold SEDs. The dashed line denotes the position of {Ve/Va) = 0.5 on the graph, values above which imply a concentration 
of galaxies at higher redshifts and below which imply a concentration of galaxies at lower redshifts. The grey shaded region shows the 
area enclosed by the Icr error. 



We also performed a banded {Ve/Va) analysis on our 
sample of AzTEC sources (excluding the AzTEC6 counter- 
part) the results of which are shown in figure [T] but our 
sample is too small to find any clear evidence of a redshift 
cutoff. 



5 CONCLUSIONS 

We give new Subaru, CFHT and IRAC photometry for a 
number of sources in the AzTEC / COSMOS survey with 
accurate coordinates from SMA imaging. We have estimated 
photometric redshifts for four of the seven galaxies in the 
sample. We find a median redshift of Zmean ~ 2.57 and 
a maximum of Zmax = 4.50. Of the sources in the com - 
bined 850 nm surv e ys presented inlChataman et al. I ll2005h 



(l2008h . consisting of ~200 sources, only one is at a redshift 
greater than our two highest redshift sources. This in ad- 
dition to the fact that we are unable to detect two of our 
sources in the optical bands down to very faint magnitudes 
may indicate that 1.1 mm surveys are more efficient at de- 
tecting very high-redshift sources than 850 fim surveys. 

Re-investigating the space densi ty evolution of a 
sample of 38 GOODS -N sources (|iPope et al. I l2006l : 
IWaU. Pope fc Scott |[2008l ') with more realistic SEDs we find 
a redshift cutoff at 2; ~ 1 if we assume a 'hot' SED and 
marginal evidence for a cutoff at z ~ 2 if we assume a 'cold' 
SED (in reasonable agreement with Wall et al.). Similar to 



Dmea oou nm surv e ys presentea miunapman et ai. 1 iizuuai. 
iPope et al. I (|2006l ). iDve et al. I (|2008l ) and IClements et aTl 



IWaU. Pope fc Scott I l|2008l) we also found evidence for two 
differently evolving sub-populations of SMGs, separated in 
luminosity, with high luminosity sources showing a less neg- 
ative evolution. 
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Figure 7. The distribution of the values of (Ve/Va) with z/ou, for the AzTEC sample. The panels on the left hand side uses the hot SED 
and the panels on the right hand side uses the cold SED. The dashed line denotes the position of (Ve/Va) = 0.5 on the graph, values 
above which imply a concentration of galaxies at higher redshifts and below which imply a concentration of galaxies at lower redshifts. 
The grey shaded region shows the area enclosed by the la error. 



We performed a similar test on the AzTEC sources but 
were unable to draw any reliable conclusions as the sample 
is too small. The GOODS-N sample is also relatively small, 
and therefore any evidence for redshift cutoffs and differently 
evolving sub-populations must be treated with caution. In 
order to harden our conclusions in general we require larger 
surveys with accurate redshifts. We would also need surveys 
taken over larger areas of sky in order to take into account 
the effects of cosmic variance. Future, larger surveys (e.g. 
with Herschel, SCUBA2) therefore will enable us to more 
robustly determine the nature of the number density evolu- 
tion of SMGs in the Universe. 
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